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Abstract. The advent of 8m-class telescopes has made galaxies at 1 ^ 
z <J 4 relatively easy to detect and study. This is a brief and incomplete 
review of some of the recent results to emerge from surveys at these red- 
shifts. After describing different strategies for finding galaxies at z ;> 1, 
Q ■ an d the differences (and similarities) in the resulting galaxy samples, I 

summarize what is known about the spatial clustering of star-forming 
galaxies at z ;> 1. Optically selected galaxies are the main focus of this 
review, but in the final section I discuss the connection between optical 
j> ■ and sub-mm samples, and argue that the majority of the 850/Uin back- 

ed \ ground may have been produced by known optically selected populations 

*0 ■ at high redshift. Among the new results presented are the dust-corrected 

luminosity function of Lyman-break galaxies at z ~ 3, the estimated con- 
tribution to the 850/xm background from optically selected galaxies at 
Q\ \ 1 ^ z ^ 5, revised estimates of the spatial clustering strength of Lyman- 

0^ ■ break galaxies at z ~ 3, and an estimate of the clustering strength of 

star-forming galaxies at z — 1 derived from a new spectroscopic sample 
Oh 1 of ~ 800 galaxies with (z) = 1.0, o~ z = 0.1. 
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Introduction 



An ambitious goal in cosmology is to understand how the universe evolved from 
its presumed beginning in the Big Bang to the familiar collection of stars and 
galaxies that we observe around us today. The last decade has seen tremendous 
progress in understanding the large role that gravitational instability almost 
certainly played. Although we still do not have complete analytic understand- 
ing, reasonable analytic approximations for the growth of gravitationally driven 
perturbations are now known, and sophisticated N-body simulators and simu- 
lations are freely available for obtaining more precise or detailed information. 
Remarkable progress has also been made in observationally constraining the ini- 
tial conditions that are required as input to the simulations or approximations. 
The available data appear largely consistent with the idea that primordial fluc- 
tuations were Gaussian (e.g. Bromley & Tegmark 1999) with a power-spectrum 
similar to that of an adiabatic ACDM model over ^> 5 orders of magnitude in 
spatial scale (e.g. numerous recent results from observations of the cosmic mi- 
crowave background, reviewed most recently by Scott 2000; White, Efstathiou, 
& Frenk 1993; Croft et al. 1999). 

But gravitational instability is only half — the easy half! — of the story. It 
alone cannot tell us how or when the stars that populate the universe today were 
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formed. Presumably stars began to form within overdensities in the matter dis- 
tribution as these overdensities slowly evolved from small ripples in the initial 
conditions into the large collapsed objects of today, but modeling this has proved 
tremendously difficult. We cannot easily model the formation of a single star 
(see Abel's contribution to these proceedings), let alone the ten billion stars in a 
typical galaxy. Even the most sophisticated theoretical treatments of galaxy for- 
mation rely on simplified "recipes" for associating the formation of stars with the 
gravitationally driven growth of perturbations in the underlying matter distribu- 
tion. The adopted recipes for star formation, although physically plausible, are 
by far the most uncertain component in theoretical treatments of galaxy forma- 
tion. We will need to check them through observations of star-forming galaxies 
at high-redshift before we can be confident that our understanding of galaxy 
formation is reasonably correct. These observations, and their implications, are 
the subject of my talk. 



2. Finding Galaxies at z ^ 1 

In the past 5 years several techniques have been shown effective for finding 
galaxies at z ^ 1. I don't have space to list them all; a partial list would include 
deep optical magnitude limited surveys (e.g. Cohen's contribution to these pro- 
ceedings), narrow band surveys (e.g. Hu, Cowie, & McMahon 1998), targeted 
surveys around known AGN at z ^ 1 (e.g. Hall & Green 1998, Djorgovski et 
al. 1999), 850/im surveys (e.g. Ivison et al. 2000), and color-selected surveys 
(e.g. Steidel et al. 1999, Adelberger et al. 2000). Different selection techniques 
have different advantages and are optimized for answering different questions. 
Color-selected surveys, which detect numerous galaxies over large and (hope- 
fully) representative volumes, are especially well suited for studying large scale 
structure at high redshift. They will be the main focus of this review. 

In color-selected surveys, spectra are obtained only for objects with broad- 
band colors indicating that they are likely to lie at a given redshift. The left panel 
of Figure 1 illustrates why galaxies at certain redshifts have distinctive broad- 
band colors. The right panel shows spectroscopic redshifts for galaxies satisfying 
various simple color selection criteria, demonstrating that color selection is a 
reasonably effective way of finding galaxies at a range of redshifts z ^ 1. At 
z ~ 1 these galaxies were selected by exploiting the Balmer-break (Adelberger 
et al. 2000), at z ~ 3 and z ~ 4 by exploiting the Lyman-break (Steidel et al. 
1999), and at z ~ 2.2 by an approach similar to the "UV drop in" technique 
described by Roukema in these proceedings. The data in Figure 1 represent only 
our own efforts; many more galaxies at similar (and higher) redshifts have been 
found by other groups with a variety of techniques. 

Although different strategies for finding galaxies at z ;> 1 result in samples 
weighted towards different types of objects, there is nevertheless significant over- 
lap between the galaxy populations that are found. The left panel of Figure 2, 
showing the Ly-a equivalent width distribution of color-selected Lyman-break 
galaxies at z ~ 3, illustrates the point. About 20% of Lyman-break galaxies 
have equivalent widths large enough to be detected in standard narrow-band 
searches for high redshift galaxies. At fixed continuum luminosity, narrow-band 
searches detect only the fraction of galaxies with the largest equivalent widths, 
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Figure 1. Left: Examples of color selection techniques. The shaded 
regions are spectra, redshifted to z = 1 and z = 3.2, of a model galaxy 
that has been forming stars at a constant rate for lGyr. Absorption 
from interstellar and intergalactic hydrogen has been applied to the 
spectrum at z = 3.2. Superimposed are the transmissivities of the 
U n GlZi filters. The Balmer and Lyman breaks give galaxies at these 
redshifts distinctive colors, allowing large numbers of them to be lo- 
cated in deep images. These and other spectral features can of course 
be used to find galaxies at other redshifts as well. Right: The redshift 
histogram of all galaxies in the color-selected samples of Steidel et al. 
(1999) and Adelberger et al. (2000). Different colors correspond to 
selection criteria aimed at different redshifts. 
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and at fixed equivalent width, color-selected surveys detect only the fraction 
of galaxies with the brightest continua; but the galaxies detected with these 
techniques appear to belong to the same underlying population. Similarly, al- 
though the z ~ 1 Balmer-break selection criteria of Adelberger et al. (2000) 
are designed to select optically bright star-forming galaxies at z ~ 1, a sub- 
stantial fraction of these galaxies (the limited available data suggests 1 in ~ 20, 
or ~ 0.2 per square arcmin to K a ~ 20) have the red optical-to-infrared colors 
1Z — K s ;> 5.5 that are often thought to be characteristic of extremely dusty 
or old galaxies at this redshift. Many of the same galaxies will therefore be 
found both by surveys for old or dusty galaxies at z ~ 1 that exploit their ex- 
pected large optical-to-infrared colors and by surveys for star-forming galaxies 
at z ~ 1 that exploit the Balmer break. Finally, there is even some overlap 
between far-UV selected samples and far-IR selected samples of galaxies at high 
redshift, though these two selection strategies might have been expected a pri- 
ori to find completely different populations of objects. For example, the two 
/^(850/im) > 5m Jy sources robustly identified with star- forming galaxies (as 
opposed to AGN) at z £ 2, SMMJ14011+0252 at z = 2.565 (Ivison et al. 2000) 
and West-MMDll at z = 2.979 (Chapman et al. 2000), have the relatively blue 
far-UV colors observed in optically selected galaxies at similar redshifts; they 
are typical, aside from their unusually bright far-UV luminosities, of the kind of 
galaxies found in optical surveys. 

The relationship between sub-mm selected and UV-selected high-redshift 
populations can be partially understood with plots like the right panel of Fig- 
ure 2, which shows the inferred distribution of dust opacities among optically 
selected galaxies at z ~ 3 (Adelberger et al. 2000). These dust opacities were 
estimated with a relationship between Ai6oo an d far-UV spectral slope that is 
obeyed by starburst galaxies in the local universe (e.g. Meurer, Heckman, & 
Calzetti 1999). It is not known if high-redshift galaxies obey this relationship; 
see §4 below. The majority of galaxies at this redshift appear to have mid- 
dling dust opacities and are therefore far easier to detect in the optical than at 
850/um, but some galaxies, even in optical surveys, are so dusty that they would 
have been easier to detect with sub-mm rather than optical imaging. Although 
relatively dust-free galaxies appear to dominate high-redshift populations by 
number, it is unclear if they dominate by star-formation rate: dusty galaxies 
tend to have much larger star-formation rates and this compensates, to some 
unknown but probably large extent, for their smaller numbers. I will discuss 
this further in 54. 



3. Spatial Clustering in High-Redshift Samples 

The large samples of star-forming galaxies at z ^ 1 produced by color-selected 
surveys allow one to begin to try to fit star-forming galaxies into the larger 
context of structure formation in the universe. The most obvious way to make 
a connection between the observed galaxies and perturbation in the underly- 
ing distribution of matter is to attempt to estimate the masses associated with 
individual galaxies by observing their velocity dispersions. Unfortunately this 
approach is surprisingly difficult. To begin with, it is hard to measure velocity 
widths for high-redshift galaxies. The [Oil], H/3, and [OIII] nebular emission 
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Figure 2. Left: The equivalent width distribution of Lyman-break 
galaxies at z ~ 3 (Steidel et al. 2000). Roughly 20% of LBGs have 
rest equivalent widths W\ ^> 20A large enough to be detected in stan- 
dard narrow band surveys. Right: The distribution of dust opacities in 
Lyman-break galaxies at z ~ 3, estimated from their far-UV spectral 
slopes j3 (fx oc \P) as described in the text. The solid bars show the 
best estimate; dotted and dashed lines indicate the size of systematic 
uncertainties. Galaxies with ^4i6oo < have broad-band colors contam- 
inated by strong Ly-a emission. Galaxies with Aiqqq ^ 5 will be too red 
to satisfy the Lyman-break selection criteria, but most known galaxies 
at z ~ 3 are much less obscured. Although relatively unobscured galax- 
ies appear to dominate the universe at z ~ 3 by number, they may not 
dominate by star formation; see §4. A galaxy with / v (850/im) = 2m Jy, 
barely detectable by SCUBA, will be too faint in the optical to be 
included in our Lyman-break survey if it has ^4i600 k; 4.5; a galaxy 
with TZ = 26, barely detectable in the deepest Lyman-break survey 
fields, will be too faint at 850/iin to be detected by SCUBA if it has 
-4i600 ^ 4.5; see Adelberger & Steidel (2000). 
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lines of galaxies at z ~ 3, for example, are redshifted into the bright sky of the 
near-IR, and as a result perhaps only 2-3 velocity widths can be measured per 
night even with an 8m-class telescope. A more serious problem is interpreting 
the velocity widths that have been measured. The limited available data suggest 
that most Lyman-break galaxies at z ~ 3 have nebular line widths corresponding 
to a 1D ~ 70-80 km/s, for example (Pettini et al. 1998, 2000). These velocity 
dispersions are far smaller than the circular velocities that were expected for 
Lyman-break galaxies on a number of other grounds (e.g. Baugh et al. 1999; 
Mo, Mao, & White 1999) and this has led some to suggest that Lyman-break 
galaxies are low mass "satellite galaxies" undergoing mergers (e.g. Somerville's 
and Kolatt's contributions to these proceedings). But because the baryons in 
these galaxies presumably cooled and collapsed farther than the dark matter 
before stars began to form, their nebular line widths are expected to be signifi- 
cantly smaller than the full circular velocity of the dark matter potential. The 
exact size of the difference is not easy to calculate. Analytic attempts at the 
calculation (e.g. Mo, Mao, k, White 1999) rely on a large number of simplifying 
assumptions, but the real situation may not be so simple: Lehnert & Heckman 
(1996) and Kobulnicky & Gebhardt (2000) have presented evidence for a com- 
plicated relationship between nebular line widths and circular velocities in the 
local universe among late-type and starburst galaxies that are presumably the 
closest analogs to detected high-redshift galaxies. 

The spatial distribution of high redshift galaxies provides an alternate way 
of making a connection between star-forming galaxies and perturbations in the 
underlying distribution of dark matter. For a given cosmogony the spatial distri- 
bution of matter at any redshift is straightforward to calculate with simulations 
or analytic approximations. Once large numbers of star- forming galaxies have 
been detected near a single redshift we can therefore ask, for example, what 
kinds of collapsed objects in the expected distribution of mass at that redshift 
have the same spatial distribution as the observed galaxies. In this way we can 
attempt to place star-forming galaxies in the larger context of structure forma- 
tion. Observational constraints on the spatial clustering of galaxies at z ^ 1 will 
be the main subject of this section, starting at z ~ 1 and moving later to z ~ 3. 

Attempts to measure clustering strength for galaxies at z ~ 1 have been 
carried out by Carlberg et al., Le Fevre et al., and Cohen et al.; their results 
are reviewed in Cohen's contribution to these proceedings. Here I will focus 
instead on previously unpublished results from a survey of star-forming galaxies 
at z ~ 1.0 (Adelberger et al. 2000). This sample consists of several thousand 
photometric candidates with 1Z <^ 25.5; to date redshifts have been obtained for 
~ 800 of them. The mean redshift of the spectroscopically observed candidates 
is (z) ~ 1.0 and the standard deviation is a z ~ 0.10. We currently have uniform 
spatial spectroscopic sampling in four 9' and one 6.5' square fields. Roughly 100 
redshifts have been obtained in each (Figure 3). 

In the four 9' fields, the variance of galaxy counts in cubes of comoving side 
length / ~ 6ft, -1 Mpc (Qm = 0.3, f^A = 0.7), estimated as described in Adel- 
berger et al. 1998 from the data in Figure 3, is (Tg al = 1.2 ±0.2. For a power-law 

spatial correlation function of the form £(r) = (r/ro) -1 ' 7 — which is consistent 
with the angular clustering of these galaxies — this variance corresponds to a 
comoving correlation length of ~ 3/i _1 Mpc, or a variance of galaxy counts in 
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Figure 3. Redshift histograms in the 5 most complete fields of the 
Balmer-break survey at z ~ 1. The field B20902 is ~ 6.5' x 6.5'; 
the others are ~ 9' x 9'. The shaded curve in the background is our 
estimated selection function, determined from the full spectroscopic 
sample. Only about a quarter of the photometric candidates were 
observed spectroscopically in each field, but the spatial sampling of 
our spectroscopy is close to uniform in these fields. 
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spheres with radius 8h 1 Mpc of cr| gal ~ 0.3. The 8h 1 Mpc variance is similar 
to the expected variance of mass at z — 1 in spheres of the same size, estimated 
by evolving back to z ~ 1 with linear theory the value of erf determined at z ~ 
from the abundance of galaxy clusters (e.g. Eke, Cole, & Frenk 1996). Balmer- 
break galaxies are evidently fairly unbiased tracers of mass fluctuations at z ~ 1 
(for the Qm = 0.3, f^A = 0.7 cosmology assumed throughout). Evolving their 
clustering strength forward to z ~ with the linear prescription of Tegmark & 
Peebles (1998) suggests that these galaxies are likely the progenitors of galaxies 
with us ~ 1 in the local universe, i.e. relatively normal galaxies. This result 
could perhaps have been anticipated from the comoving abundance of Balmer- 
break galaxies, which, at ~ 0.02/i 3 Mpc~ 3 to 1Z ~ 25.5, is similar to that of L* 
galaxies in the local universe. 

The increase in clustering strength from galaxies in our sample at z ~ 1 
to galaxies at z ~ therefore appears to be relatively easy to understand; it is 
almost exactly what one might have expected gravitational instability to have 
produced acting on a population of formed objects. Remarkably the same is not 
true at higher redshifts. Rather than decreasing further, at redshifts z > 1 the 
observed clustering strength of detected star-forming galaxies begins to rise. 

Hints that star-forming galaxies at z ;> 2 might be strongly clustered were 
first provided by targeted surveys of small and carefully selected volumes, often 
around known AGN (e.g. Giavalisco, Steidel, & Szalay 1994; Le Fevre et al. 
1996; Francis et al. 1997; see also later work by Campos et al. 1999 and 
Djorgovski et al. 1999). Further evidence came subsequently from the color- 
selected survey of Lyman-break galaxies at z ~ 3 (e.g. Steidel et al. 1998). 
Figure 4 shows the projected correlation function Lo p (r p ) (e.g. Davis & Peebles 
1983) of galaxies in this sample. The implied correlation length, neglecting 
systematic errors, is 3.8±0.3/i _1 Mpc comoving (Om = 0.3, Q\ = 0.7). A similar 
estimate of the correlation length follows from the relative variance of Lyman- 
break galaxy counts in cubes of comoving side-length 11.4/i _1 Mpc, <7g al = 0.75± 
0.25 (Adelberger et al. 2000; this value supersedes our previous estimate, which 
was based on a smaller data set). 

Gravitational instability acting on a population of galaxies that had ro ~ 
4/i -1 comoving Mpc at z ~ 3 would not produce a population with ro ~ 3/i _1 
comoving Mpc (similar to the observed ro of Balmer-break galaxies) at z ~ 1 or 
a population with ro ~ 5/i _1 comoving Mpc (similar to normal local galaxies) at 
z ~ 0. This can be shown in a crude way by first assuming that the correlation 
function of Lyman-break galaxies selected at z ~ 3 would maintain a constant 
slope of 7 ~ 1.7 at lower redshifts, and then using the linear approximation 
of Tegmark & Peebles (1998) to evolve the observed clustering of Lyman-break 
galaxies at z ~ 3 to lower redshifts. Limited space does not allow a more careful 
analysis or the consideration of cosmological models besides Q.m = 0.3, $7 a = 0.7, 
although both would affect our conclusions somewhat. In this simplified analysis, 
we would expect former Lyman-break galaxies to have ro ~ 5/i _1 comoving Mpc 
at z ~ 1 and ro ~ 9/i -1 comoving Mpc at z ~ 0. 

These correlation lengths are significantly larger than those of Balmer-break 
galaxies at z ~ 1 and of normal galaxies at z — 0, implying that Lyman-break 
galaxies are unlikely to be the progenitors of either population. What are they 
instead? Why is their spatial clustering so much stronger than we might naively 
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Figure 4. The projected correlation function of Lyman-break galax- 
ies at z ~ 3, from Adelberger et al. (2000). The quoted confidence 
limits on ro and 7 do not include systematic uncertainties, which are 
likely at least as large as the random uncertainties. 



have expected? A clue is provided by their number density, ~ 4 x 10~ 3 /i 3 Mpc" 3 
(fijv/ = 0.3, Q\ = 0.7) to 7Z = 25.5, which is about 5 times lower than the number 
density of normal galaxies in the local universe or of Balmer-break galaxies with 
1Z ^ 25.5 in our z — 1 sample. Perhaps the relatively rare Lyman-break galaxies 
are not progenitors of typical galaxies at z ^ 1 but instead are special in some 
way. One way in which they are special is that they are the UV-brightest 
galaxies (and therefore presumably the most rapidly star-forming galaxies) at 
z ~ 3. Semi-analytic calculations (e.g. Baugh et al. 1998, Kauffmann et al. 
1999) suggest that the most rapidly star-forming galaxies at high redshift will 
reside within the most massive collapsed objects, rather than within typical 
collapsed objects, and so perhaps we can understand the clustering of Lyman- 
break galaxies by trying to associate them with massive collapsed objects at 
z ~ 3 instead of with galaxy populations at lower redshift. 

In a classic paper, Kaiser (1984) showed that in hierarchical models the 
most massive collapsed objects at any redshift are more strongly clustered than 
the distribution of matter as a whole. A formalism for estimating the clustering 
strength of collapsed objects as a function of mass was subsequently developed 
by many authors; see (e.g.) Mo & White (1996). Remarkably the observed 
clustering of Lyman-break galaxies is indistinguishable (as far as we can tell — 
see Wechsler's contribution to these proceedings) from the predicted clustering 
of the most massive collapsed objects at z ~ 3 down to a similar abundance (e.g. 
Adelberger et al. 1998). This result suggests that there may indeed be a simple 
relationship between mass and star-formation rate in high redshift galaxies, as 
many semi-analytic models predicted (e.g. Baugh et al. 1998). 
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If the mass of a galaxy plays a dominant role in determining its star for- 
mation rate, then we might expect the star-formation rate distribution of z ~ 3 
galaxies to be related in a simple way to the distribution of masses of collapsed 
objects at z ~ 3. As a crude guess at the star- formation rate associated with 
a collapsed "halo," we can take the mass cooling rate in the halo for large 
masses and the mass cooling rate times a number proportional to v 2 c for small 
masses where supernova feedback is important (e.g. White & Frenk 1991). In 
this approximation, for cooling dominated by Bremsstrahlung, we would expect 
SFR « M 5 / 6 for large M and SFR Z M 3 / 2 for small M. Figure 5 shows the 
slopes of our simplistic "theoretical" SFR distribution in these two limits. At 
each abundance the slope of the mass function was estimated with the Press- 
Schechter (1974) approximation for an Qm = 0.3, Q\ = 0.7, h = 0.65, a§ = 0.9, 
r = 0.2 cosmogony; changing the values of any of these parameters would change 
the mass function slope somewhat. Also shown in Figure 5 is the observed "dust- 
corrected" luminosity function of these galaxies, estimated using the (3-Aiqqq 
correlation of Meurer et al. (1999) as described in Adelberger Sz Steidel (2000). 
The observed slopes agree reasonably well with our naive expectations, but the 
star-formation rates at any abundance are unexpectedly high, many times larger 
than the expected cooling rate for halos of similar abundance. (Standard for- 
mulae can be used to derive star-formation rates from the far-UV luminosities 
shown in Figure 5; see, e.g., Madau, Pozzetti, & Dickinson 1998.) This has been 
taken as evidence that the star formation in Lyman-break galaxies is fueled not 
by quiescent cooling but instead by the rapid cooling that would accompany 
a merger of two smaller galaxies (e.g. Somerville's and Kolatt's contributions 
to these proceedings) . Many uncertain steps lie behind the conclusion that the 
star formation rates in Lyman-break galaxies far exceed their quiescent cooling 
rates, however. The star-formation rates of Lyman-break galaxies could be con- 
siderably lower than is usually deduced, for example, if we are wrong about the 
shape of the IMF, about the magnitude of the required dust corrections, or even 
about the value of various cosmological parameters, while the cooling rates in 
Lyman-break galaxies could be significantly higher than usual estimates if we 
are wrong about the baryon fraction, the metallicity, or the spatial distribution 
of the gas or dark matter in these galaxies. It would be tremendously interest- 
ing if the star formation rates of Lyman-break galaxies really were far higher 
than their quiescent cooling rates, but I do not think this has been conclusively 
shown. 

In any case I will assume for now that Lyman-break galaxies are strongly 
clustered because they reside within rare and massive collapsed objects at z ~ 3, 
and return to the question of how they might be related to the Balmer-break 
galaxies observed at z — 1. The calculation above showed that z ~ 3 Lyman- 
break galaxies with 1Z ^ 25.5 cannot be the progenitors of z ~ 1 Balmer-break 
galaxies with 1Z ^ 25.5, and this was hardly surprising since the abundance of 
the Lyman-break galaxies is so much lower. But suppose we had significantly 
deeper U n GlZ photometry, so that we could detect fainter Lyman-break galaxies 
and reach an abundance similar to that of Balmer-break galaxies at z ~ 1. Could 
this deeper population of Lyman-break galaxies evolve into a population like the 
Balmer-break galaxies by z ~ 1? If halo mass and star formation rate are related 
in the simple way described above, then the deep Lyman-break population would 
be somewhat less strongly clustered that the current 1Z <J 25.5 population, help- 
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Figure 5. The observed and dust-corrected luminosity functions of 
Lyman-break galaxies at z ~ 3. At this redshift the 1Z filter samples ~ 
1700A rest, and so the (dust-corrected) 1Z luminosity should be closely 
correlated with star-formation rate. Error bars do not include some 
systematic uncertainties; these uncertainties are especially severe for 
"^•corrected 24 — see Adelberger & Steidel (2000) before using this plot. 
The solid lines show the expected slope of a naive "theoretical" star- 
formation rate (SFR) distribution for SFR oc M (left) and SFR cx 
M i.5 ( r ight). See text. 



ing to remove the inconsistency between the observed tq of Balmer-break galaxies 
and the expected ro of Lyman-break galaxy descendants. A deep Lyman-break 
population, with abundance ~ 15 times that of the 1Z <^ 25.5 sample, has in fact 
been detected in the HDF, and its correlation length ro (which cannot be mea- 
sured very accurately because of the HDF's small size) appears smaller than that 
of the brighter ground-based population (Giavalisco et al. 2000); but it looks like 
this effect is not strong enough to remove the inconsistency in clustering strength 
for Lyman-break galaxy descendants and Balmer-break galaxies at z ~ 1. We 
are left with the result that galaxies selected by the Balmer-break technique at 
z ~ 1 are probably not (for the most part) the descendants of those detected 
with the Lyman-break technique at z ~ 3. Because Balmer-break galaxies ap- 
pear to be representative of typical star- forming galaxies at z ~ 1, the simplest 
interpretation is that Lyman-break galaxies have largely stopped forming stars 
by z ~ 1 — are they perhaps instead passively evolving into the elliptical galaxies 
observed at lower redshifts? 

In principle this sort of argument could provide a very stringent constraint 
on the star- formation rates of Lyman-break galaxy descendants at z ~ 1, since 
(for Qm = 0.3, S7a = 0.7) galaxies with as little as ~ O.2/i~ 2 M /yr of star 
formation would have bright enough UV continua (in the absence of dust obscu- 
ration) to be included in our Balmer-break sample. But I am still not convinced 
that the argument is completely robust; firmly establishing that former z ~ 3 
Lyman-break galaxies are not significantly forming stars at z ~ 1 will require 
a much better understanding than is currently available of the clustering of 
fainter and more numerous Lyman-break galaxies at z ~ 3 and of what kinds of 
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star-forming objects might not satisfy our "Balmer-break" selection criteria at 
z ~ 1. 



4. Dust 

A full observational understanding of star formation at high redshift can only 
be achieved if we are able to detect most of the star formation in representative 
portions of the high redshift universe. But how can we be sure that our surveys 
are not missing a large fraction of the star formation at high redshift? The sad 
answer is that we can't. There are too many ways that star formation could be 
hidden from our surveys for us ever to be sure that we have detected most of it. 
Surveys cannot detect the star formation that occurs in objects below their flux 
and surface brightness limits, for example, and they cannot directly detect the 
formation of the low mass stars that (at least in the local universe) dominate 
the total stellar mass. At best we can aim for a reasonably complete census of 
the formation of massive stars that occurs in objects above our flux and surface 
brightness limits. If these limits are deep enough, and if the high-redshift IMF 
is similar enough in all environments to what we have assumed, then this sort of 
sample can serve as an acceptable proxy for a true census of all star formation 
at high redshift. 

What is the best way to produce a reasonably complete survey of massive 
star formation? Massive stars emit most of their luminosity in the UV, and 
so naively we might choose a deep UV-selected survey. It has become clear in 
recent years, however, that most of the UV photons emitted by stars in rapidly 
star- forming galaxies are promptly absorbed by dust, and as a result most of 
the luminosity produced by massive stars tends to emerge from these galaxies 
in the far-IR where dust radiates. This is true in the local universe for a broad 
range of rapidly star-forming galaxies, from the famous class of "Ultra Luminous 
Infrared Galaxies" (ULIRGs, galaxies with Lfir ^ 10 12 L Q ) to the much fainter 
UV-selected starbursts contained in the IUE Atlas (e.g. Meurer et al. 1999); 
and the recent detection of a large extragalactic far-IR background (e.g. Fixsen 
et al. 1998) suggests that it is likely to have been true at high redshifts as well. 

Two implications follow. The first is that far-IR luminosities probably pro- 
vide a better measurement of rapidly star- forming galaxies' star- formation rates 
than do UV luminosities. The second is that even very rapidly star-forming 
galaxies may not be detected in a UV selected survey if they are sufficiently 
dusty. Far-IR/sub-mm selected surveys should therefore in principle provide 
a much better census of massive star formation at any redshift than UV se- 
lected surveys. 850/um surveys are likely to provide an especially good census 
at 1 <; z <J 5, because favorable K-corrections make a galaxy that is forming 
stars at a given rate appear almost equally bright at 850/um for any redshift in 
this range (see Hughes' contribution to these proceedings), and consequently a 
flux-limited sample at 850/um is nearly equivalent to a star-formation limited 
sample at 1 <J z <J 5. 

This is exactly the kind of sample that detailed attempts to understand star 
formation at high redshift require. So why then was most of this review devoted 
to UV-selected high redshift samples rather than the apparently superior 850/im 
samples? The reason is that 850^m observations are comparatively difficult. The 
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deepest 850/mi image taken with SCUBA (the current state-of-the-art sub-mm 
bolometer array) reached a depth of 2mJy in a ~ 4 square arcminute region of 
the sky, for example, and five sources were detected (Hughes et al. 1998). In 
contrast a modern instrument on a 4m-class optical telescope can easily obtain 
photometry to a depth of <J 0.1/xJy over a ~ 500 square arcminute region, 
detecting thousands of galaxies at z ^ 1. Even though optical surveys do not 
select star-forming galaxies in the optimal way, most of what we know about 
galaxies at high redshift comes (and will continue to come for many years) from 
these surveys. The question — probably the most important question for those 
interested in high-redshift galaxy formation — is whether the large and detailed 
view of the high redshift universe provided by UV selected surveys is reasonably 
complete. 

This is equivalent to asking whether the galaxies responsible for producing 
the 850/um extragalactic background are bright enough in the rest-frame UV to 
be included in current optical surveys. If they are not, then optically selected 
surveys will not be able to teach us much of value about high-redshift star 
formation despite the wealth of information they contain. 

The straightforward way to constrain the UV luminosities of the objects 
that produce the 850/um background is to observe the UV luminosities of known 
850/um sources. This is difficult in practice, however, because relatively few 
850/um sources have been detected and each has a significant positional uncer- 
tainty due to SCUBA's large diffraction disk. Often several optical sources lie 
within a sub-mm error box, and a great deal of effort is required to determine 
which one is the true optical counterpart (e.g. Ivison et al. 2000). Moreover only 
about 30% of the 850/um background can be resolved into discrete sources with 
current technology, and so optical observations of detected 850/um sources can- 
not conclusively tell us about the optical luminosities of the galaxies responsible 
for producing most of the 850/mi background. 

An alternate approach, taken by Adelberger & Steidel (2000), is to estimate 
the contribution to the 850/xm background from known optically selected popu- 
lations at high redshift and compare this expected background to the observed 
background to see if there is significant shortfall. 

Meurer et al. (1999) have shown that the far-IR luminosities of UV-selected 
starbursts in the local universe can be estimated to within a factor of ~ 2 from 
the starbursts' UV luminosities and spectral slopes j3. Their /3/far-IR relation- 
ship was the foundation of our calculation. It is not known if high-redshift 
galaxies obey this relationship. Chapman et al. (2000) presented evidence that 
the 850//m fluxes of z ~ 3 Lyman-break galaxies may be somewhat lower than 
the relationship would suggest, but the evidence is very marginal when the uncer- 
tainties in Lyman-break galaxies' predicted 850/um fluxes are taken into proper 
account. Adelberger & Steidel (2000) have checked in a number of other ways 
whether z ^ 1 star-forming galaxies obey the relationship. Figure 6, showing 
the predicted and observed 15//m fluxes of Balmer-break galaxies in the z ~ 1 
sample described above, is an example. The predicted fluxes in Figure 6 as- 
sume that these galaxies will follow both the /3/far-IR correlation of Meurer et 
al. (1999) and the correlation between 6-9 ;um luminosity and far-IR luminos- 
ity observed in the (star- forming) ULIRG sample of Genzel et al. (1998) and 
Rigopoulou et al. (2000). The 15/xm data are from the ISO LW3 observations of 
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Figure 6. The predicted and observed 15/im fluxes of Balmer-break 
galaxies at z ~ 1. Fluxes were predicted from these galaxies' UV prop- 
erties using the /3/far-IR correlation of Meurer et al. (1999) and the 
6-9/mi/far-IR correlation observed among rapidly star- forming galax- 
ies in the local universe. The majority of these galaxies are undetected 
at a 3cr limit of 150/xJy; uncertainties on the predicted fluxes of these 
galaxies have been omitted for clarity, and the upper limits to their ob- 
served fluxes are indicated with downward pointing triangles. The six 
detections are indicated by points with error bars in both x and y; the 
one with Si5 !pr edicted *C Si5,observed has a significant separation between 
its optical and mid-IR centroids and may be a misidentification. 



Flores et al. (1999). A large fraction of objects with predicted fluxes below the 
detection limit were not detected, and a large fraction of objects with predicted 
fluxes above the detection limit were detected. This plot therefore provides some 
support for the notion that Balmer-break galaxies at z ~ 1 obey the /3/far-IR 
relationship. 

If we assume that all optically selected galaxies at z ^ 1 obey this relation- 
ship, then we can make crude estimate of their total contribution to the 850/mi 
background. In this calculation I will assume further that the comoving star 
formation density in optically selected populations is constant for 1 < z < 5 
(e.g. Steidel et al. 1999), that the dust SEDs of optically selected galaxies at 
1 < z < 5 are similar to those of starbursts and ULIRGs in the local universe, 
and that the (unknown) luminosity and (3 distributions of optically selected 
galaxies at 1 < z < 5 are similar to those measured for Lyman-break galaxies at 
z ~ 3. (See Adelberger & Steidel 2000 for a more complete discussion.) Under 
these assumptions optically selected populations at 1 < z < 5 would be expected 
to produce 850/um number counts and background that are surprisingly close to 
the observations (Figure 7). 

Although the overall agreement is good, within the substantial uncertain- 
ties, there appear to be significant differences at the brightest 850/um fluxes. 
Optically selected galaxies cannot easily (it seems) account for the large num- 
ber of observed sources with /^(850/xm) ^ 6mJy, and indeed Barger, Cowie, 
& Richards (1999) have shown that these sources tend to have extremely faint 
optical counterparts. It is possible that these bright sources are associated with 



15 




Figure 7. Top panel: The contribution to the 850^m background 
from UV-selected populations, calculated as described in the text. Con- 
straints on the bright end from Barger, Cowie, & Sanders (1999) are 
shown as points with uncertainties. Our best guess at these galaxies' 
background contribution (the area of the shaded curve) agrees to within 
10% of the measured value from Fixsen et al. (1998), but the excel- 
lent agreement is likely only a coincidence; the systematic uncertainties 
in this calculation are very large. The dashed lines illustrate the sys- 
tematic uncertainty due to possible errors in our derived (3 distribution 
alone. Bottom panel: Observed and predicted 850/mi cumulative num- 
ber counts. The observed number counts, from Blain et al. 1999, are 
denoted by points with error bars. Our best guess prediction is the 
solid line; dashed lines are as above. 



AGN, rather than the star-forming galaxies we have included in our calculation, 
but in any case it appears that the bulk of the 850^m background could have 
been produced by known optically selected populations at high redshift. This 
claim rests on a large number of assumptions that could easily be wrong, but I 
am willing to bet $100 that when ALMA finally resolves the ~ mm background 
we will discover that most of it is produced by galaxy populations already de- 
tected and studied in the rest-frame UV. Any takers? 
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